Dissipation scales of kinetic helicities in turbulence 
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A systematic study of the influence of the viscous effect on both the spectra and the nonlinear 
fluxes of conserved as well as non conserved quantities in Navier-Stokes turbulence is proposed. This 
analysis is used to estimate the helicity dissipation scale which is shown to coincide with the energy 
dissipation scale. However, it is shown using the decomposition of helicity into eigen modes of the 
curl operator, that viscous effects have to be taken into account for wave vector smaller than the 
Kolomogorov wave number in the evolution of these eigen components of the helicity. 

PACS numbers: 47.65.-d, 52.65Kj, 91.25Cw 



I. INTRODUCTION 

In two recent papers it was suggested that dissipation 
of kinetic helicity occurs at a scale kjj 1 larger than the 
Kolmogorov scale k^ 1 . This was justified on dimensional 
grounds \M as well as using a GOY shell model of tur- 
bulence [2(. In contrast using a different shell model of 
turbulence based on helical wave decomposition, both 
scales were found to be equal kn = ks Q- I n addition, 
direct numerical simulations, also presented in 0], seem 
to confirm the latter result though, as noted by the au- 
thors, the computational limitations prevent to have a 
Reynolds number sufficiently large to really discriminate 
between both scenarii. 

The purposes of the present work are to investigate 
further the possible existence of a specific helicity dissi- 
pation scale and to understand why two shell models do 
exhibit different helicity behaviours while their energy 
spectrum are very much similar. Part of this apparent 
contradiction comes from the very definition of the dissi- 
pation scale. Indeed, in the Kolmogorov theory, there is 
no ambiguity. The scale at which the energy dissipation 
terms are no longer negligible when compared to the non- 
linear fluxes of energy corresponds to the scale at which 
the energy spectrum departs from the Kolmogorov power 
law. This scale marks the end of the cascade process as 
well as the beginning of energy spectrum fall off. 

The situation is less clear for non conserved quanti- 
ties such has the positive H + and the negative helicity 
H ~ defined respectively as the helicity carried on by the 
eigenvectors of the curl operator with positive and nega- 
tive eigenvalues. Generally, for non conserved quantities 
Q, we propose to refer to the dissipation scale as the scale 
after which the dissipative term dominates the dynamics 
so that the spectrum of Q falls off. Such a scale might 
very well differ from the scale, referred hereafter as the 
viscous scale, at which the dissipative terms start to play 



a role in the dynamics of Q. Indeed, for a non-conserved 
quantity, the nonlinear term might very well compen- 
sate for the increase of dissipation in part of the high 
wave number range after the viscous scale and prevent 
the spectrum to fall off even if dissipation is active. In 
general, the viscous scale should be smaller than the dis- 
sipation scale. However, for conserved quantities, both 
the viscous and the dissipation scales coincide. 

A general discussion on the determination of dissipa- 
tion scale is presented in Section II for conserved as well 
as non conserved quantities. The specific case of the two 
conserved quantities in three-dimensional turbulence, the 
energy and the helicity, is discussed in Section III. The 
positive and negative helicities, which are not conserved 
quantities, are discussed in Section IV. Shell models de- 
scribing the high Reynolds number behaviour of turbu- 
lence are discussed in Section V. Both models used in 
references (@) and (||) are introduced and analyzed nu- 
merically in Section IV. It is shown very clearly that the 
dissipation scales for the helicity and the energy coincide 
and are given by the Kolmorogov length scale. More- 
over, the dissipation scale for the positive and negative 
helicities also corresponds to the energy dissipation scale. 
However, the analysis of their fluxes allows to identify 
very clearly a viscous scale for both H + and H~ that is 
smaller than the energy dissipation scales. 



II. DISSIPATION SCALES IN TURBULENT 
SYSTEMS WITH CASCADES 

Before discussing the specific problem of energy or he- 
licity dissipation scale, we consider a general quadratic 
quantity Q that is not necessarily conserved by the non- 
linearities of the Navier-Stokes equation: 



= / d 3 r a(r)6(r) = / d 3 k o(k)6(k) 



c.c. . 



(1) 



Here, a(r) and 6(r) are two fields and a(k) and 6(k) arc 
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assumed to be statistically isotropic. In that case, it is 
convenient to introduce the spectrum Q(p), so that: 



Q = J dp Q( P ). 



(2) 



The parts of this quantity that are represented by modes 
such that |k| < k and |k| > k are denoted respectively 
by Q<(k) and Q>{k): 



Q < {k)= / d 3 k a(k)6(k)* + c.c, 

J|k|<K 

Q > (k) = [ d 3 k o(k)6(k)* + c.c. , 

J|k|>K 

Q = Q < (k) + Q>(k) Vk. 

Their evolution is given by: 

d t Q<( K ) = s Q -U<( K )-d<( K ) 
d t Q>(n) = -U>(n)-d>^) 



(3) 

(4) 
(5) 

(6a) 
(6b) 



where sq is the source of Q here injected by a forcing 
process in the largest scales of the system (kp ) so that 
kp < k. In that case, the source term is independent of 
k. The nonlinearity contributions to the evolution of Q < 
and Q > arc noted respectively Uq(k) and IIq(k). They 
correspond to fluxes respectively outward and inward the 
sphere of radius n if Q is a conserved quantity. The 
dissipation of Q in the modes |k| < k (resp. |k| > k) is 
noted dg(n) (resp. dq(n)). 

In the following, the dissipative processes are assumed 
to be represented by viscous type terms, so that: 

d<(K) Q = 2v f dpp 2 Q(p) (7a) 
Jo 

/>oo 

d>{K) Q = 2v dpp 2 Q{p). (7b) 



If the system undergoes a cascading process that trans- 
fers Q from the forcing scales to small scales, the non- 
linear transfer at scale k should be characterised by a typ- 
ical time scale that will be denoted Tq(k). On the other 
hand, dissipation processes should also be characterized 
by a time scales Tq 1ss (k). In the case of viscous type dis- 
sipation, Tq 1ss (k) = 1/(vk 2 ). The comparison of these 
characteristic time scales can be used to estimate the end 
of the cascade range (usually referred to as the inertial 
range as long as kinetic energy is concerned). Indeed, 
in the range dominated by the non-linear interactions, 
t q( k ) < t q ss ( k ) since non-linear interactions should be 
faster than dissipative processes. On the contrary, in the 
dissipation range Tq(k) > Tq ss (k). An estimate of the 
dissipation scale kg is thus: 



Hiss ( iD 



(8) 



Of course, in order to predict kg, it is necessary to guess 
the expression for Tq(k). For instance, if a scaling law 



can be assumed, Tq(k) = Aq k aQ , the dissipation wave 
number is given by: 



k D 



1/(2-cq) 



(9) 



Another typical length scale can be introduced via 
Eq. (pja)) and under the assumption that a stationary state 
can be reached: 



dpp 2 Q(p). 



(10) 



This expression can be used to obtain an estimate of 
the viscous scale kq at which the viscous term becomes 
important when compared to sq by assuming that the 
spectrum Q(p) follows a power law: Q(n) = Bq k~@ q : 



dpp 2 B Q p 



which leads to 



"■■r> 



1/(3-/3q) 



(11) 



(12) 



For a conserved quantity, the spectrum has to fall off for 
K > k V Q otherwise the dissipation would exceed the in- 
jection rate and consequently the nonlinear term must 
vanish. It is thus expected that kq = kq. However, 
for a non conserved quantity, the dissipation may exceed 
the injection rate since the nonlinear term does not nec- 
essarily vanish. Thus, nothing prevents the spectrum to 
remain Q(k) = Bq k~@ q after kq and kq may be smaller 
than kq. 



III. ENERGY AND HELICITY DISSIPATION 
SCALES 

We first consider the cascade of energy. The total en- 
ergy injection rate is then usually noted se = e and the 
Kolmogorov energy spectrum can be derived, 



E(k) = C E e 2 / 3 A;" 5 / 3 



(13) 



in the inertial range. The estimate for the dissipation 
wavenumber based on the equality between the charac- 
teristic time scales requires an expression for t e 1 (k). Var- 
ious proposals can be found in the literature, but all yield 
the same scaling since they are built with only k and e, 
assuming the viscosity does not influence the non-linear 
characteristic time: 

T E \ K ) (X K^/V 1 / 3 (14) 

which means Ae x r _1 /' J and <\ E = 2/3. Consequently, 
expression ^ yields 



oc 



c-1/3 



3/4 



1/4 



(15) 
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Similarly, the Kolmogorov spectrum implies (3e = 5/3 
and B E oc e 2 / 3 , and the expression (fl"2"|) yields the same 
estimate: 

/ e x 3/4 , £ s 1/4 

We now consider the helicity cascade. Both studies 
presented in [l[ and Q make the assumption that the 
characteristic time of non-linear transfer of energy and 
helicity are the same: Tjj(k) = (k). Since, both en- 
ergy and helicity arc dissipated by linear viscous pro- 
cesses, their dissipation characteristic time are obviously 
identical t^ ss (k) = rf ss (/«) = l/(vk 2 ). In that case, 
the dissipation wavenumber for energy and helicity ob- 
tained by comparing the non-linear transfer time to the 
dissipation time must coincide: 

1/4 



kjj oc kj<} oc 



(17) 



Also, the equality of the non-linear transfer time is also 
known to imply the following helicity spectrum : 

H{k) = C H Se- 1/3 k-^ 3 , (18) 

where Cr is a dimensionless constant and 8 is the helicity 
injection rate. In that case, the formula (fT"2)) with sjj = S 
and Bh = St^ 1 / 3 leads to the same expression 



k H oc 



3/4 



e \ 1/4 



(19) 



Hence, both approaches yields the same result and tend 
to confirm the equality between the helicity and the en- 
ergy dissipation scales. However, although the equality of 
both dissipation scales is so obvious, the analysis becomes 
a bit more involved when using the helical decomposition 
of the energy and helicity spectra. 



IV. HELICAL DECOMPOSITION OF SPECTRA 

Following the approach presented in [l[, the Fourier 
modes of both the velocity and the vorticity are expanded 
using a basis of polarised helical waves defined by 
ikxh^ ±kh ± 0: 

u(k) = u+(k)h+ +u-(k)h _ , (20) 

u>(k) = fcu+(k)h+ - ku-(k)h~. (21) 

The energy and helicity carried on by the mode u(k) 
respectively become 

u(k)-u*(k)/2 = (| U +(k)| 2 + | U -(k)| 2 )/2, (22) 
u(k)-w*(k)/2 = fc(|u+(k)| 2 - |u-(k)| 2 )/2. (23) 

Isotropy is again assumed and both the energy E(k) 
and the helicity H{k) spectra are considered to be func- 
tions of k = |k|. Introducing the spectral densities of 
energy and helicity for the helical modes (±) yields: 

E{k) = E + (k)+E-(k), (24) 
H{k) = H + {k) + H~(k) = k[E+{k) - E-(k)} .(25) 



Their equations of evolution have exactly the struc- 
ture ©. Moreover, all these quantities are dissipated 
through viscous effect and their linear dissipation time 
scale is again r dlss (K) = \/(vk 2 ). Guessing their non- 
linear characteristic time is however much more difficult. 
Indeed, non-linear transfers can transform E <+ (k) not 
only in E >+ (k) but also in E < ~(n) and E > ~(n). More- 
over, E + and E~ are not separately conserved by the 
non-linear terms. Hence, invoking the equality of char- 
acteristic time scales to estimate the dissipation scales of 
these quantities is not really possible. 

It is however quite easy to estimate their spectra 
from (fl"3l) and (fTSl): 



E +(k) = ^ £ 2 / 3 fc" 5 / 3 + ^( ( y/ e i/3)j fc -8/3 ) (26a ) 

E -(k) = ^-e 2 ' 3 ^' 3 - ^(5/e^ 3 )k^ 3 , (26b) 

which are the equations (9) and (10) of [l|. As a conse- 
quence, the leading order in k must be given by 

£±(fc) = ^ £ 2/3 fc -5/3 5 H ± {k) = ± ^ £ 2/3 fc -2/3_ ( 27) 



By construction, the range of validity of (]26[) and (|27|) 
is the same as that of the scaling laws (fT3| and (|18[) of 
E(k) and H(k). It is therefore bounded by = kj^. 

On the other hand, formula ([T2"j) yields an estimate 
of the scale from which on the dissipative term must be 
considered in the evolution of £ ,± and _ff ± . It leads to 



k 



E± 



1/4 



and k H ± 



c2/3 



3/7 



6 3 



1/7 



(28) 



(29) 



As noted by Ditlevsen and Giuliani [lj, k u H± < k E ± = 
kg. Indeed, the helicity injection rate is at most /cp e, so 
that 



1/7 



k D 



3/7 



K E ! 



(30) 



and since kp < kg in the turbulent regime, k u H± < kg. 
However, there is no reason to identify the helicity dissi- 
pation scale as k H± . Clearly, the spectrum of H + {k) and 
H~{k) can not deviate from the scaling k~ 2 l 3 and fall off 
in the range k H± < k < k E . Indeed, if these quantities 
decay faster that k~ 2 / 3 after k H '± , then the quantities 
E ± (k) = ±H ± (k)/k will decay faster than fc" 5 / 3 and 
k v H± would be identified as the end of the inertial range 
which is known to actually extend down to k E . Consid- 
ering the Eq. (fT0|) applied to Q = H*- allows to better 
understand the meaning of k H ± : 



n±< = S ± t ^C E ve 2 l 3 k 7 ' 3 . 



(31) 



For low values of k <C k v H± , the fluxes are constant and 
equal to 6 . However, for k ^> k H ±, the dissipation of 
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H ± is stronger than the injection rate 5^ and the non- 
linear flux has to scale like fc 7 / 3 . Ditlevsen and Giuliani 
referred to this scale the dissipation scale of Q ± . How- 
ever, as argued above, this does not correspond to the 
end of the helicity spectrum. 



V. HELICAL SHELL MODEL ANALYSIS 

Shell models are built to describe the exchange of phys- 
ically relevant quantities between the various scales of a 
turbulent flow. The Fourier space is divided into a set 
of shells which are logaritmically spaced. A field (like 
velocity for instance) is represented by very few (1 or 2) 
complex variables in each shell. These models allow to 
investigate turbulence properties at a much lower numer- 
ical cost than direct numerical simulation. Actually, the 
resolutions achievable in DNS are still too limited to dis- 
tinguish clearly k u H± from Q. On the other hand, as 
will be done below, Reynolds numbers - defined at scale 
kp as Re = e 1 l 5 k F i ^' i v~ 1 - as large as 10 7 can be reached 
with a shell model. 

Helical shell models were developed in and are 
based on the helical decomposition of Fourier modes [|[ . 
Two dynamical variables are used per shell to repre- 
sent both the helical components of the velocity field. 
As pointed out in [(|, such helical models can be re- 
trieved from the helical triadic systems of the Navier- 
Stokes equations in helical basis. Four simple models 
can be expressed in a single formula: 



d t u ri 



+ /, 



± 

n ! 



(32a) 



with 



Wt = ik n [{ Sl \-s 2 \ 2 )u^ + \u± 



+ 2 



± SI ±S2S1 



+ (s 2 A — A *) u. 
+ (A- 2 - SiA" 1 ) ut\< ai 



n-l "n+1 

±s 

n-l 



(32b) 



where each model is obtained for one particular choice of 
(si, S2) with si, S2 = ±1. In (|32() the parameter A is the 
logarithmic shell spacing and the wave number is defined 
as k n — /uq \ n . 

In the absence of forcing and viscosity u, the shell 
model (f32|) conserves total energy E and helicity H @ . 



JV 



jV 



E — ^i?„, H — H n , 



(33) 



71=1 



where N is the number of shells in the model. The energy 
E n and helicity H n in shell n are defined as 



-E^ 
Hi 



E- 



Et 



1 



,±|2 



±|2 



(34) 
(35) 



Within the model, the fluxes of energy and helicity are 
defined as 



n<(n) 



r< 1 



n£<(n) 



■ni<(n), 



n^(n) = n+<(n) + H H < 
with the following explicit expressions: 



(36) 
(37) 



n|<(n) 



n|<(n) = 



71 

= -J2W±ui* + ce, (38) 

m— 1 

- [dtjz\(±k m )\uf n A 

\ m=l J NL 

n 

= TJ2 k mW±ut* + cc, (39) 



pled sets of variables appear, namely: [u{ 
and 



u 2 1 u 3 1 ■ 



m—l 

where ^lg(k) is the flux (due to the non- linear term) of 
the quantity Q leaving the region of wave numbers lower 
than k and Hq < (/s) is the flux leaving either the '+' or 
the '— ' variables of wave numbers lower than k. 

The GOY model used in Q corresponds to (si,S2) = 
(— , +) in the helical picture (|32|). In this case, two uncou- 

. ). In the original version of the GOY 
model, only one of these sets is considered. Hence, in 
each shell n, the helicity is evaluated alternatively by 
H% or H~, depending whether n is odd or even. The 
cancellation of the leading terms in equation ()25[) with 
the scaling (|26a[) and (|26b[) does not occur. Therefore 
H(k) can not be straightforwardly obtained with a GOY 
model. The fluxes presented in are hence closer to 
n#(K) than to n^(re) although, stricto sensu they are 
neither of them. 

On the other hand, the developments proposed in Q 
were illustrated by the SABRA version of the model cor- 
responding to (si, S2) = (+, — ) in which all variables are 
coupled. Both H£ and H~ are available within each 
shell n and so is the total helicity H n . In the following 
section, the work of Q is pursued and the energy and 
helicity spectra and fluxes are investigated. 



VI. NUMERICAL RESULTS 

The computation of the averaged helicity spectra, 
which is the difference of its two helical components and 
requires the cancelling of the leading terms demands very 
fine time stepping. Furthermore, very long simulations 
are required in order to obtain enough statistics. This is 
probably the reason why helicity spectra have not been 
reported so far [D-Q, with the notable exception of 0. 
Very long and accurate integration of the shell model (|32|) 
with (si, s 2 ) = (+, — ) have been performed. In these sim- 
ulations, the forcing is concentrated on one single shell 
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FIG. 1: Helical case: e = e + — 1, S~ — e~ — 0. Energy and 
helicity plots are respectively represented on the left and right 
columns versus logfc. The positive and negative helical modes 
are denoted by o and •, the sum of both modes by +. The 
spectra (resp. fluxes) are represented in the top (resp. down) 
row. 



(the fourth) and provides constant energy and helicity 
injection rates. The rate of energy injection within the 
'±' variables is denoted e . The rate of helicity injec- 
tion is therefore 5 ± = ± kp£ ± ■ The phases are randomly 
chosen at each time step. 

In figures Q] and [2] the results are presented for respec- 
tively a helical and a non helical case. The parameters 
are v = 10" 7 and A = (1 + \/5)/2. The shell are la- 
belled from —2 to 37 with k n — A™. The total number 
of shells is thus N = 40 and the forcing is concentrated 
in the third shell so that fej? = 1. For the helical case 
e + = e = 1, implying S = S + = 1 and e~ = —5~ jkp = 0. 
For the non helical case e + = e~ = 1/2, implying e = 1, 
S + = —5~ = 1/2 and 5 = 0. In each figure, the left and 
right columns correspond respectively to energies and he- 
licitics. 

The spectra are plotted in log- log frames (upper row). 
Energies E(k n ) and E ± (k n ) scale in fc^ 2 ^ 3 correspond- 
ing to power spectral densities in fc -5 / 3 in agreement 
with (Ti~3l and (|2"T|) . Helicities H ± (k n ) scale in fc„ cor- 
responding to power spectral densities in fc -2 / 3 in agree- 
ment with (f2~T)) . In the helical case, the total helicity 
H(k n ) scales in k n 2 ^ 3 corresponding to a power spectral 
density in fc~ 5 / 3 in agreement with (fTS)) , In the non he- 
lical case H(k n ) is the sum of two opposite quantities 
H ± (k) and has no clear scaling. Compared to H ± {k n ) it 
can be considered as negligible, in agreement with (TT5)) 
taking 8 = 0. Note that all spectra manifestly extend up 
to the Kolmogorov scale k^ ~ 10 s 

The nonlinear fluxes are plotted in log-log frames 
(lower row). For the helical case, the total energy flux as 
well as the energy flux of E + are constant and dominated 
by e = 1 up to the Kolmogorov scale. On the contrary 
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FIG. 2: Same as figure\T\for the non helical case e = e + = 
1/2, 5 = 0. 



the flux of E~ has no component corresponding to the 
injection since e~ = so that its spectrum is dominated 
for low k by the viscous term and is proportional to fc 4 / 3 . 
The viscous scale k v H± is clearly identified on the helic- 
ity flux for H + . For k < k v H± , the flux is constant and 
dominated by 8 + while for k < k v H± the injection is sub- 
leading and the flux scales like k 7 ^ 3 . Remarkably, the 
viscous scale k < k v H± is also very clearly observed even 
in the non-helical case. 



VII. CONCLUSION 

The present study has allowed to identify two different 
length scales related to to the dissipation of a quadratic 
quantity Q in Navier-Stokes turbulence. The first one 
is the traditional dissipation scale that marks the end of 
the power law in the spectrum of Q due to the domi- 
nant effect of the viscosity. The second scale, referred 
to as the viscous scale, corresponds to the beginning of 
the range in which viscous effect have to be taken into 
account. Clearly, for the kinetic energy, the viscous and 
the dissipation scales coincide. However, for non con- 
served quantities, such as the positive and negative part 
of the helicity, these two scales arc different. Although 
the viscous scale cannot be measured from the spectra, 
it is easily identified from the nonlinear fluxes. This has 
been shown using shell models. 

This approach reconcile the analysis of [l], HJ and 0] . 
Strictly speaking, the scale k v H± cannot be interpreted as 
the dissipation scale for helicity. Both direct shell model 
integration and helical components analysis show that 
the helicity cascade develops down to the Kolmogorov 
scale. However, this scale is indeed relevant in the anal- 
ysis of the nonlinear flux of helicity and plays a role even 
when the flow is globally non helical. 
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Beyond the issue of hclicity dissipation scale which is 
now clarified, this study stresses how much caution is re- 
quired when studying the effect of helicity on turbulence 
dynamics with a GOY model Q ■ Other models like the 
one used here or those presented in @ or Q are highly 
preferable. 
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